Covalency effects on the magnetism of EuRh2P2 
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In experiments, the ternary Eu pnictide EuRh2P2 shows an unusual coexistence of a non-integral 
Eu valence of about 2.2 and a rather high Neel temperature of 50 K. In this paper, we present 
a model which explains the non-integral Eu valence via covalent bonding of the Eu 4f-orbitals to 
P2 molecular orbitals. In contrast to intermediate valence models where the hybridization with 
delocalized conduction band electrons is known to suppress magnetic ordering temperatures to at 
most a few Kelvin, covalent hybridization to the localized P2 orbitals avoids this suppression. Using 
perturbation theory we calculate the valence, the high temperature susceptibility, the Eu single- 
ion anisotropy and the superexchange couplings of nearest and next-nearest neighbouring Eu ions. 
The model predicts a tetragonal anisotropy of the Curie constants. We suggest an experimental 
investigation of this anisotropy using single crystals. From experimental values of the valence and 
the two Curie constants, the three free parameters of our model can be determined. 

PACS numbers: 71.10.-w, 71.27.-fa, 75.30.Ds, 75.50. Ee 



I. INTRODUCTION 

Non-integral valence, intermediate valence, and 
magnetic order 

A general question of interest is how a non-integral 
valence of localized ions in a solid influences the pos- 
sibility of magnetic order. Concerning this subject, an 
earlier paper reported the anomalous valence state of 
Eu in EuRh2P2'i EuRh2P2 was characterized to show 
intermediate-valent and probably also covalent proper- 
ties. A coexistence of a non-integral Eu valence of 2.2 
and antiferromagnetic order up to 50 K was reported. 

Simply given, a non-integral valence means that the 
mean total occupation number of the ionic electronic lev- 
els is non-integral. This might have several physical rea- 
sons. In particular, intermediate valence is the hybridiza- 
tion of localized ionic states with the strongly delocalized 
conduction band. It is experimentally evident and the- 
oretically well-understood that there is a strong compe- 
tition between intermediate valence and magnetic order. 
Intermediate valence is known to suppress magnetic or- 
dering temperatures typically to at most a few K. As an 
experimental example, TmSe is intermediate-valent and 
has a Neel temperature below 3.5 Ki^ 

A detailed theoretical approach to intermediate va- 
lence is the extended s-f modeli^ Using this model a simi- 
larity to the Kondo effect is shown: grossly spoken, infor- 
mation about the electron spin is washed out by the de- 
localization of the conduction band states. The extended 
s-f model allows one to take a ferromagnetic exchange be- 
tween ionic f and conducting s states into account, which 
is typical for intermediate-valent Eu ions. Hence, it is 
more realistic than an Anderson model, which always 
implies an antiferromagnetic s-f couplingi^ 

In agreement with experiment, the extended s-f model 
also explains why the magnetic ordering temperature of 
an intermediate-valent system is enhanced drastically by 
mechanical or chemical high pressure. In addition, a 



given intermediate valence is pressure-sensitive itself be- 
cause the energy levels of the localized states can be 
pressed towards the conduction band. These aspects 
are not too important for the present paper because we 
will concentrate neither on intermediate valence nor on 
anomalous pressure. 

If a non-integral valence is caused by covalence, no 
strongly delocalized states are involved in the underlying 
hybridization. Hence, a characteristically different ap- 
proach will be required to understand such systems. 

The interpretation of EuRh2P2 is insufficient so far. On 
the one hand, if EuRh2P2 was a system in which a Eu in- 
termediate valence of 2.2 and magnetic order coexist up 
to 50 K without any restrictions, this would mean a small 
sensation. Not only was such a behaviour experimentally 
unusual but also theoretically not understood, because in 
this case the ordering temperature was an order of mag- 
nitude higher than for typically realistic parameter sets 
of the extended s-f model. '^ Even for extreme choices of 
the model parameters (in case of a Eu intermediate va- 
lence of 2.2), an upper boundary of about 15 K for the 
ordering temperature is estimated. A new concept of 
intermediate valence would have to be found to explain 
the experiment. On the other hand, no detailed model is 
available to investigate the counter-perspective, the in- 
fluence of the covalence on the magnetism of EuRh2P2. 

In the following, we will motivate the use of a covalent 
scenario instead of intermediate valence as the starting- 
point for the investigation of EuRh2P2. This is in con- 
trast to former publicationsi:^:^^ and also refers to un- 
published experimental materia L^'^'^^i^^ 



The covalence of EuRh2P2 

We take into account the following experi- 
ments on EuRh2P2: measurements of the magnetic 
susceptibilityjii^ of the crystal structure r'iii Lm x-ray 
absorption,* and Mofibauer spectroscopyii^ 




FIG. 1: Crystal structure of EuRh2P2 



The properties of EuRh2P2 can be unraveled by study- 
ing the influences of temperature and doping, especially 
with As. Depending on these two parameters, in exper- 
iments a structural phase transition of first order is ob- 
served whereas the lattice type, a body-centered tetrag- 
onal ThCr2Si2 structure (figure[T]), does not change. 

The structural phases a and f3 differ by a strongly 
anisotropic jump of the lattice parameters and by the 
electronic bonding conditions. This leads in the f3 phase 
to a non- integral Eu valence of 2.2. The lattice cells show 
a slight elongation in the square plane but a considerable 
compression along the tetragonal axis^iiiii compared to 
the a phase, which contains divalent Eu ions. In both 
phases, the valence is homogeneous, i.e., all Eu ions are 
equivalent. 

In the P modification two nearest neighbouring P 
atoms form a single molecular bondliii — oriented along 
the tetragonal axis — which does not appear in the a 
phase. That corresponds to a charge transfer from the 
phosphorus to the conduction band, which leaves holes 
in P2 molecular states. These are available for occupa- 
tion by one fluctuating Eu electron each. In this way, in 
the P phase there is a covalent hybridization between the 
magnetic Eu ions and P2 molecules. The Rh ions do not 
participate in the magnetic properties of EuRh2P2. 

Both phases a and /3 exhibit a magnetic phase transi- 



tion from para- to antiferromagnetism at a Neel temper- 
ature Tn (phase diagram: figure 2). With decreasing As 
doping Tn drops somewhat at the a-/3 phase boundary 
but does not change its order of magnitude. According to 
the extended s-f model this effect of the valence transition 
on the ordering temperature is too small to be consistent 
with intermediate valence. 

The valence measurements on EuRh2P2, in particular 
Lni x-ray absorption and MoBbauer spectroscopy, can 
not distinguish between a non-integral valence of cova- 
lent and intermediate-valent origin. The reduced mag- 
netic moment of the P phase compared to the divalent 
Eu moment of the a phase does not clarify the nature of 
the non-integral valence either. Furthermore, the valence 
cannot be derived from the magnetic Eu moment alone 
because in the non intermediate-valent case a non-trivial 
contribution by the P2 states must be considered. This 
quantity is unknown. 

The covalent bonding scenario between Eu and P2 
states persued in the present paper was derived in ref- 
erence [Il| from bonding lengths and is supported by fur- 
ther aspects. In particular, the pressure dependences^ 
of the valence, which we express as the deviation AW 
from the divalent state, leads to the conclusion that the 
intermediate-valent part of the total Eu valence amounts 
to at most 



AW^Iiv « 0.05. 



(1) 



This is the precision with which the Eu valence is deter- 
mined by MoBbauer spectroscopy. ^° Reference 10 states 
no intermediate valence within this precision because in 
MoBbauer spectroscopy the Eu valence does not change 
over the range between ambient pressure and 5 GPa. 
This is the expectation in absence of a conduction band 
hybridization. 

As the results of Lni x-ray absorptioii^' may be influ- 
enced by final-state effects, we assume an error bar for 
the non-integral Eu valence of the /3 phase including in- 
termediate valence as well as covalence and estimate the 
valence shift as: 



AM/ = 0.20 ± 0.05. 



(2) 



We will not carry out an error calculation in closed 
form but will select several discrete values for parame- 
ters of the model, which will be constituted neglecting 
AH^Iiv. Furthermore, we ignore the weak temperature 
dependence of the valence which is considerably less sig- 
nificant in the /? phase than the precision of 0.05.^° 

The mechanism and the geometric structure of mag- 
netic ordering in EuRh2P2 are almost completely un- 
known. Experimentally, antiferromagnetism below 
Tn has been concluded, as well as the existence 
of — unspecified — ferromagnetic couplings because of an 
anomalous paramagnetic Curie temperature.— In addi- 
tion, the importance of a phosphorus-mediated superex- 
change between the Eu ions has been shown qualitatively 
via the sensitivity of the MoBbauer magnetic hyperfine 
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FIG. 2: Phase diagram of EuRh2(Asj;Pi_3;)2.® The mid-grey 
sector corresponds to a coexistence due to hysteresis. 



turbatively from divalent ionic ground states we view 
the crystal as built by inactive two-dimensional metal- 
lic Rh planes and by quasi-two-dimensional insulating 
EuP2 planes. The covalent hybridization in these planes 
is much more important than the one between different 
EuP2 planes. This is concluded from the Eu-P distances, 
which in the (3 phase near the phase boundary are: 3.88 
A between the planes and 3.10 A in the plane. 

As the starting-point we view the properties of an 
arbitrary single Eu ion in interaction with the four 
neighbouring P2 molecular ions which form a tetrago- 
nal cage around it. This interaction describes the high- 
temperature paramagnetism of the crystal to leading or- 
der. We begin with three unknown model parameters 
and reduce their number finally to one after calculating 
two quantities for which experimental values are avail- 
able: the valence and the paramagnetic susceptibility (as 
an avarage over the three spatial directions) . We also cal- 
culate the single-ion anisotropics of the Eu ions, which 
are experimentally unknown. 

Using the same model parameters, since thermal ex- 
pansion is very small in the temperature range between 
and about 400 K)i we calculate the superexchange pa- 
rameters between nearest and next-nearest neighbouring 
Eu ions. 



field (and Tn) to pressure,— which is not observed in the 
reference system EuRh2As2. 

We conclude that the covalence of EuRh2P2 has to 
be significant because the intermediate-valent part is 
small and the reduced magnetic moment, the existence of 
the superexchange, and the bonding scenario have been 
shown to be mutually consistent with the covalence. 

A possibility other than covalence or intermediate va- 
lence which can cause a measured non-integral valence 
should be excluded in EuRh2P2: In Eu(Pdo.7Auo,3)2 Si2 
similar Lni x-ray absorption and M6f5bauer results in 
spite of a divalent magnetic Eu moment are found at 
ambient pressure and are explained finally by an anoma- 
lous spatial extension of the 4f shell, whereas the sys- 
tem becomes intermediate-valent under high pressure. ^^ 
In contrast to EuRh2P2, in Eu(Pdo,7Auo,3)2 Si2 no co- 
valent bonding partners are available for the Eu ions. 

Under pressure beyond 5 GPa, intermediate valence in 
EuRh2P2 acquires, at least qualitatively, the same rele- 
vance as the covalence but stays less important than the 
latter at ambient pressureJ^ The corresponding charac- 
teristic decrease of Tn under increasing pressure is un- 
derstood by the extended s-f model. 

Band structure calculations of EuRh2P2 have only been 
performed for the intcgral-valent a (high temperature) 
phasci^ 

The magnetic ordering of EuRh2P2 is caused both by 
super- and indirect exchange. The latter aspect follows 
because the reference system EuRh2As2 exhibits almost 
the same T^ (48 K) but neither superexchangei^ nor an 
anomaly of the paramagnetic Curie temperature^S, 

In order to treat the covalence effects of EuRh2P2 per- 



il. MODEL 



The a phase of EuRh2P2 has got the formal valencesii: 
Eu2+(Rh2+)2 (P3")2 . (3) 

full 3p shell 

This corresponds to the divalent ground state of a Eu 



4f^ S^^, L^O, J 



(4) 



In the covalent (3 phase there is a fluctuation between the 
formal valences)^ 



Eu2+(Rh^ 



Eu3+(Rh+)^ (P2)'-(5) 



1 hole per 
P core 



1 hole per 
molecule ion 



The hybridization of the Eu ions is de- 
scribed perturbatively — starting from the divalent 
configuration — via a quantum mechanical admixture of 
trivalent states: 



4f^ 5" = 3, L^3, J = 0...6. 



(6) 



Hund's rule correlations are taken fully into account. 
We use a Lande approximation for the energies Ej of the 
Eu'^^ ground state J — and the low excitations J > 0: 



Ej = AE[l + XJ{J+l)]. 



(7) 



Ej^o is shifted by an unknown charge transfer energy 
AE with respect to the divalent ground state, which is 



the zero point in our calculation {E2+ = 0). The in- 
traionic spin-orbit part of the levels Ej is given by 



X 



e 



A2AE 



^ = 7960K 



(8) 



where I5M) is a state of the Eu^+ configuration (M = 
7/2. ..7/2), \JJz) is a state of the Eu'^"'" configuration 



(J = 0...6,J^ 
element. 



J... J), and f . , is a reduced matrix 



WjJ 



which was fitted to optical measurements^ by a least 
squares fit. ^ is the spin-orbit coupling parameter {Hgo = 
^LS is the Hamiltonian of the spin-orbit coupling for a 
single Eu ion). 

After an electron hops from the Eu ion to a phosphorus 
ion, it occupies an anti-bonding molecular orbital (MO) 
state which is considered as a linear combination of two 
3pz atomic orbitals due to the MO method and which is 
odd under a reflection with respect to the a^ {xy) plane. 
The bonding 3p^ MOs are always filled. 

As the perturbation term which expresses the covalent 
hybridization we introduce a hopping operator between 
a Eu ion and four neighbouring P2 molecular ions: 



v^J2 [tofL + i-^)''t,{.fUa+fL) 



H'' 



H. 



(9) 



fc^l,2,3,4 



p and / denote the annihilators corresponding to the 
single-particle states. Because of symmetries (time re- 
versal of the crystal Hamiltonian, reflection with respect 
to xy and xz plane) contributions due to hopping am- 
plitudes t„i for magnetic quantum numbers 771 ^ 0, ±2 
are excluded and we have t-2 — i2- Pa relates to the 
antibonding MO of the fcth neighbour. There are three 
unknown model parameters: X, io? and ^2- The hopping 
amplitudes can be chosen real. 

We calculate the matrix elements of effective operators 
via standard perturbation theory for a degenerate system 
and to leading order in the perturbation. We use a formu- 
lation due to Takahashii^ involving a unitary transforma- 
tion r which maps the perturbed problem onto the un- 
perturbed ground state space. T is given as a power series 
in terms of the unperturbed Hamiltonian — corresponding 
here to the energies E2+ = and Ej — and the pertur- 
bation operator V . Any operator A in the Hilbert space 
of the perturbed states is treated as an effective operator 
a = r^Ar in the ground-state space. 

In the calculation, the matrix elements of the effective 
operators are expressed in terms of the matrix elements 
of the / creators (which appear in the perturbation term 
V) between the correlated many-body states of the Eu^+ 
and Eu'^+ ions. The matrix elements are evaluated via 
Clebsch-Gordan coefficients and the Wigner-Eckart the- 
orem. Explicitly, using Clebsch-Gordan coefficients the 
/ creators are transformed from the ma to the jjz basis: 



fLa = y! {im.\(j\jjz)f]. 



33^ 



(10) 



The Wigner-Eckart theorem gives: 

{W\fkVJ^) == \\f\\jj{JJ^n.\lM), (11) 



III. SINGLE-ION EFFECTS 

We calculate the matrix elements of the effective opera- 
tors of the valence and the magnetization due to hopping 
of second order in V , which is the leading order of the 
hybridization for the single-ion effects. For this calcula- 
tion it is sufficient to apply the unitary operator T to first 
order. 

We have calculated the effective single-ion Hamilto- 
nian to second order in V . According to this calcula- 
tion, which we do not present in detail, the octets of 
the unperturbed Eu^^ ions split into four Kramers dou- 
blets which, however, remain quasi-degenerate, i. e. the 
splitting of these doublets is very small compared to the 
temperature. Hence, we will use the average energy of 
the Kramers doublets for thermodynamic averaging. 



Effective valence 

The covalent admixture of Eu^+ states shifts the va- 
lence from 2 to a larger value. In order to calculate the 
model valence we use the valence operator: 



W ^'i-Pa. 



(12) 



Pq is the projector onto the unperturbed ground-state 
space J^ In the framework of our perturbation methodic 
we use the effective valence operator w = T'^WT to sec- 
ond order in V . Thermodynamic averaging due to the 
mean eigenvalues of the effective Hamiltonian of second 
order gives the mean valence: 



{W) 



■ trw. 



(13) 



Our calculation of the deviation AW = (W) — 2 from 
the divalent Eu configuration caused by the covalence 
gives the result: 



AW ^ — 



t§ 



2i| 



{AEy 



E 



1-1-2 J 



[l + J(J-f 1)X]2 



(14) 



Because of the uncertainty of the precise value of AW 
[see ^ and ([2|)] we carry out the model calculation using 
the experimental values 



AW = 0.2 and AW = 0.15. 



(15) 



Fixing the value of AW accordingly reduces the num- 
ber of unknown model parameters from three to two. For 
convenience, we define 



t 



ti 



2ti 



(16) 



as the total hopping amphtude. Because of i_2 = ^2, 
t'^ is proportional to the total hopping probability of all 
single- particle Eu states m involved in the hopping. In 
view of Eq. ([H]) , t can be considered a function of X at 
a given value of AW. 

In order to ensure that the second-order calculation 
be consistent, t/AE must be considerably lower than 1. 
This will be the case for our choices of the model pa- 
rameters. Higher orders will be suppressed by a factor of 
{t/AE)^. 



Effective Curie susceptibility 

Due to the covalent admixture of the Eu'^^ configura- 
tion the paramagnetic susceptibility of EuRh2P2 can be 
expected to become anisotropic. According to the space 
group the susceptibility tensor has tetragonal symmetry. 
In the high-temperature regime we are considering here 



the susceptibility is Curie-like. In the limit of small dif- 
ferences between the energies of the Kramers doublets, 
the paramagnetic susceptibility of the system per space 
direction i — x,y, z is given by: 



Xi: 



kT ^i--' 



^ii — 



IE 

MM' 



(|Af|m,||Af')|', (17) 



following from the standard formula for the Curie suscep- 
tibility of single ions.— The Curie constants Cu of second 
order depend on the effective magnetic moment rui = 
r^ Afi r of second order in V, where Mi = Ji + Si denotes 
the untransformed moment. Experimentally, the Curie 
susceptibility has been measured on polycrystalstii^ This 
corresponds to the spatial average C = {2Cxx + Czz)/5, 
and we define 

ACu = C-Cu (implying AC = -iAC,,). (18) 



The resulting model Curie constants are (in units of /i|): 



C 



21- A(^_ 



)^[16- 



^C..^UA)'[2A- 



+ 



1+2X 
45 

6 

1+2X 

63 



1+6X 
65 

26 
1+6X 

77 



1 + 12X 

70 
(1+12X)^ 

56 



i+i2Js: 

56 



30 

i+20Js: 

54 
(1+20X)^ 

72 
1+20X 

77 



55 


-f 


91 


1+30X 


1+42X 


11 


■- 


65 


(1+30X)^ 


(1+42X)^ 


22 




182 



i+30Js: 

143 



1+42X 



(1-|-2X)^ ' (1+6^)'' (1+12X)-' (1+30X)'' (1+42X)-' 



r 



(19) 



In this calculation we have included (in contrast to an 
intermediate-valent hybridization) non-trivial contribu- 
tions to the Curie susceptibility from the localized Pz 
orbitals: 



C^C 



CP% CP^ 



6 AW. 



(20) 



As t/AE can be expressed via X [see Eq. ([T4| ]. using 
Eqs. (|14|) and (fT9|) we can fix t and X (and consequently, 
AE) from the experimental values of AW and C. Sim- 
ilarly to the valence, we choose two experimental values 
of the Curie constants of the polycrystalline samples:- 



C = 17.2 and C = 17.6. 



(21) 



These values refer to two different samples. The latter 
value of C is inconsistent with AW^=0.2 according to the 
model. Probably, this is not an objection to the model 
but a further hint at AW < 0.2 as far as the valence 
shift is caused by the covalence. We fix three sets of 
model parameters X (or AE) and t according to Table 
m Notice that the perturbation parameter (t/AE)'^ is 



considerably lower than 1, which shows that the low- 
order calculation is sufficient. The first parameter set is 
the least realistic one because t is extremely high. 

Table U shows that for a given value of AW, the total 
hopping amplitude t is very sensitive to the value of C. 
Hence, it is important to take into account the contri- 
butions CP^ from the P2 molecules to the paramagnetic 
susceptibility, see Eq. (^0)1 . 

For convenience we define the relative hopping ampli- 
tude concerning to = single-particle states: 



to 
t ' 



(22) 



After t and X (and AE) have been fixed using the 
experimental values of AW and C, r is the only unknown 
parameter of the model. Notice that the anisotropy of 
the Curie constant, ACzz, is proportional to r^ [see Eq. 
(fT9)) ]. i. e., this anisotropy is solely caused by the covalent 
hybridization which involves m = single-particle states 
of the Eu-f shell. 




FIG. 3; Quasi two-dimensional EuP2 plane with a superex- 
change cluster of nearest neighbouring Eu ions and their com- 
mon P2 molecular neighbours. 



Our model characterizes intervals (upper bounds) for 
the anisotropy of the Curie constant. Czz/Cxx ~ 1 (> 0) 
is almost exactly proportional to r^ . Table H] lists the 
maximum anisotropy due to: 



CI 



= C, 



r~i mm r~^ 



(23) 



Single-ion anisotropy 



The effective Hamiltonian to second order in the per- 
turbation y gives rise to a single-ion anisotropy of the 
Eu ion: 



K; 



A*CT^ 



(24) 



Here and in the following we normalize every spin op- 
erator with respect to 1 as cr, — 2S,/7. The single- 
ion anisotropy parameter /i is given in Table [HI The 
anisotropy changes sign for small values of r. We will 
come back to the single-ion anisotropy in the next section 
when we discuss various effects on the magnetic ordering 
temperature, which are present in EuRh2P2. 



IV. THE COUPLING OF NEIGHBOURING EU 
SPINS 

Nearest neighbours 

The superexchange dynamics of nearest neighbouring 
Eu spins in EuRh2P2 is described to leading order by 
processes of fourth order in the EU-P2 hopping within a 
cluster as shown in Figure 3. 



In treating the intermediate states of the perturbation 
series, we use the same Lande levels as in section 2 and 
ignore Coulomb repulsion within the P2 molecules in case 
of doubly occupied P2 orbitals. In analogy to Eq. ([H]), 
the hybridization operator is given by: 



^=EK/.if + (-i)'=*2[/^^^' 



2cr 



I^] 



}py 



{k) 



H. 



(25) 
I is a site index for the Eu ions. We calculate the effective 
superexchange Hamiltonian of fourth order in V which 
scales with t^(A_E)~'^ and describes any superexchange 
in the crystal to leading order. Superexchange processes 
between nearest neighbouring Eu ions mediated by hop- 
ping paths exceeding the cluster we consider are at least 
of the order ^^(AE;)-^. 

We express the matrix elements of the effective su- 
perexchange Hamiltonian in terms of the same quanti- 
ties as in section 2. The result is a finite polynomial in 
terms of spin operators. To a good approximation the 
superexchange Hamiltonian turns out as an xxz model: 



hxxz — jx{cix<^2x + criyCr2jy) + jz ClzCr2z 



(26) 



There are additional (multilinear, see below) parts 
of the superexchange Hamiltonian, e.g., terms Cf1^a2i, 
whereas the tetragonal symmetry of the superexchange 
Hamiltonian is exact due to the crystal symmetry. The 
hopping processes we consider have got complicated se- 
lection rules. For instance, exchange processes of fourth 
order in V are possible where the z component of one 
Eu spin is changed from —3/2 to 7/2, i.e., these pro- 
cesses have the selection rule lS.Sz = 5. This is why the 
exchange Hamiltonian is not bilinear in the spin opera- 
tors but multilinear. For a bilinear exchange Hamilto- 
nian the selection rule /S.Sz = ±1 is required. However, 
the exchange Hamiltonian we obtain is bilinear to a good 
approximation. The values of the coefficients of the ne- 
glected part of the Hamiltonian depend on the choice of 
the parameter sets: less than 10"'^ ja;, 5 x lQ~^jx: and 
IG'^ix for X = 0.00181, 0.00989, and 0.0177, respec- 
tively. 

The — antiferromagnetic — coupling constants are listed 
in terms of approximate numbers in Table Hill The pre- 
cision of these numbers decreases with increasing value 
of X as the neglected parts of the Hamiltonian get more 
important. The lower the values of X and r, the more 
isotropic is the xxz coupling. The contribution, which 
is proportional to t^t^, is negligible for every parameter 
set. Figure 4 shows the coupling, which has a consider- 
able strength in each case, in terms of the spatial average 
J-(2j.+j.)/3. 



Next-nearest neighbours 

The calculation of the coupling of next-nearest neigh- 
bours (these are located within the EuP2 planes) is car- 
ried out using almost the same formalism as for nearest 
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FIG. 4: Average coupling constant of nearest neighbours. 



neighbours. The only difference is the single P2 molecule 
ion involved in this case. 

The resulting superexchange Haniiltonian for next- 
nearest neighbours is denoted by: 



Kxz = J'x{<^lx^2x + '^ly'^2!y) + fz '^lz'^2z 



(27) 



The calculated coupling constants are listed in Table 
IIVI In contrast to h^xz , now the xx part j'^ of the cou- 
pling is approximately independent of r whereas the Ising 
part j^ is lower than jz exactly by a factor of four. 



Competing effects on the magnetic ordering 
temperature 

As the value of the model parameter r, the coupling be- 
tween different EuP2 planes and the effect of the indirect 
exchange (mediated by delocalized conduction band elec- 
trons) on the intraplanar spin couplings are not known, 
we cannot present a quantitative calculation of the mag- 
netic ordering temperature. However, in the following 
we will argue why the considerable value of 50 K for the 
ordering temperature is generic in the framework of our 
model. 

There are competing effects on the magnetic ordering 
temperature of EuRh2P2. (i) There is a twofold frus- 
tration of the magnetic interactions caused by the an- 
tiferromagnetic next-nearest neighbour couplings in the 
EuP2 planes and by the coupling between different EuP2 
planes irrespective of the sign of the Heisenberg coupling 
j inter between Eu ions in different planes, see Figure 5. 
The frustration tends to decrease the magnetic ordering 
temperaturcii^ (ii) On the other hand, there may be ef- 
fects that tend to enhance the ordering temperature. Ex- 
cept for certain small values of r, one of these effects is 




FIG. 5: The frustration of the Eu lattice, caused by antifer- 
romagnetic intraplanar Heisenberg couplings j. An example 
of four antiferromagnetically ordered spins in a EuP2 plane 
is shown, (i) There is a frustration caused by the antiferro- 
magnetic next-nearest neighbour couplings j'^ and j'^ within a 
EuP2 plane, (i) A nearest neighbour belonging to a different 
plane is geometrically frustrated independently of the sign of 

J inter . 



the single-ion anisotropy (described by the parameter /x, 
see Eq. [51] and Table |TT| , which may have considerable 
strength compared to the superexchange couplings and 
may be more relevant than the xxz anisotropics of these 
couplings. The ordering temperature is enhanced drasti- 
cally, i. e., logarithmically by a single-ion anisotropy J^ 
The unknown Heisenberg exchange coupling j inter or 
an unknown dipolar interaction between different EuP2 
planes has — on the mean-field level — no effect on the or- 
dering temperature as a hypothetically given Neel order 
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in one plane would not cause a mean field on a Eu site in a 
neighbouring plane. However, the unknown anisotropics 
of the interplanar exchange and pseudodipolar couplings 
between neighbouring planes tend to stabilize the mag- 
netic order ji^ 



Ordering temperature of an isolated quasi 
two-dimensional EuP2 plane 

Though quantitative estimates of various magnetic 
couplings which may have effects on the magnetic or- 
dering temperature as discussed above are missing, we 
can present estimates for that temperature according to 
the coupling in the quasi two-dimensional EuP2 planes. 
This can serve as the starting point for a more compre- 
hensive analysis in the future, in particular including the 
interplanar couplings which are not known to date. 

Ref. [13 takes into account the isotropic Heisenberg 
couplings j of nearest-neighbouring magnetic ions on a 
two-dimensional square lattice as well as the single-ion 
anisotropics fi and estimates the magnetic ordering tem- 
perature based on these two parameters for fi << j: 



Tn = 



12j 

|ln|f| 



(28) 



This Equation gives already considerable ordering tem- 
peratures except for very narrow windows of our un- 
known model parameter r, namely Tn > 200 K except for 
r < 10-12, 0.13 < r < 0.15 and 0.21 < r < 0.22 for the 
model parameter sets with X = 0.00181, X = 0.00989 
and X = 0.0177, respectively. 

Altogether the discussion of the competing effects on 
the ordering temperature and of the intraplanar estimate 
([28|) shows that the experimentally observed ordering 
temperature of about 50 K is consistent with our cal- 
culations because in this sense our model avoids the sup- 
pression of the ordering temperature which is implied by 
intermediate-valent modelsi^ 



have argued why the experimentally observed ordering 
temperature is generic, because for instance the single- 
ion anisotropics might have considerable strength. 

TABLE I: Selected sets of model parameters and maximum 
Curie anisotropy. 



AW 


C 


X 


AE 


t 


y-i max / /^ niin 


0.2 


17.2 


0.00181 


104687 K 


24402 K 


1.018 


0.15 


17.6 


0.00989 


19163 K 


4505 K 


1.065 


0.15 


17.2 


0.0177 


11083 K 


2915 K 


1.109 



TABLE II: Single-ion anisotropy parameter. 



X 

0.00181 
0.00989 
0.0177 



14 r^ 
-1+51 r^ 
-4+87 r^ 



Measurements of the anisotropy of the Curie constants 
could fix the last free parameter r of the single-ion 
anisotropy and supcrexchange model and determine the 
model completely. Following that, measurements of the 
magnetic structure and the magnetic excitations via neu- 
tron scattering could make a description possible which 
also takes into account indirect exchange between the Eu 
ions. (The reader should be reminded that neutron scat- 
tering requires the particularly expensive isotope Eu-153 
because the standard isotope Eu-151 absorbs neutrons 
too strongly.^) In this way, there is the chance to under- 
stand the ordering mechanism in EuRh2P2 comprehen- 
sively. 
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V. SUMMARY 

We have introduced a systematic interpretation of the 
Eu valence shift and the magnetism of EuRh2P2 due 
to covalent bonding which — in contrast to a hypothet- 
ically given intermediate valence — is consistent with ex- 
periment. We have presented a model for the covalence 
which predicts upper bounds of the anisotropy of the 
Curie constants and which characterizes the strength of 
the Eu single-ion anisotropics and of the supcrexchange 
coupling between nearest and next-nearest neighbouring 
Eu ions. Though a quantitative calculation of the mag- 
netic ordering temperature has not been possible, we 
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TABLE III: Anisotropic coupling constants of nearest neigh- 
bours. 



X i. [K] 

0.00181 1091 r'' + 1091(l-r 2)2 

0.00989 133r^+ 130(1 -r2)2 



0.0177 



89r'*+ 83(1 -r2)2 



J. [K] 
1094r*+1091(l-r2)2 

136r''+ 130(l-r2)2 

94r^+ 83(1 -r2)2 



TABLE IV: Coupling constants of next-nearest neighbours. 



X 


ii [K] 


j'. [K] 




0.00181 


273 


273 r 4+273 (1 - 


-T^f 


0.00989 


33 


32r'*+ 34(1- 


-T^f 


0.0177 


22 


21r''+ 23(1- 


-r'f 
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